Diquark-Antidiquark with open charm in QCD sum rules 
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Using the QCD sum rule approach we investigate the possible four-quark structure of the recently 
observed charmed scalar mesons Dg(2308) (BELLE) and Dq' + (2405) (FOCUS) and also of the very 
narrow D+.(2317), firstly observed by BABAR. We use diquak-antidiquark currents and work to 
the order of m s in full QCD, without relying on l/m c expansion. Our results indicate that a four- 
quark structure is acceptable for the resonances observed by BELLE and BABAR: 7?g(2308) and 
D+j(2317) respectively, but not for the resonances observed by FOCUS: D°' + (2405). 
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In general, the classification of mesons containing a 
single heavy quark is interpreted with the help of heavy- 
quark symmetry, i.e., the symmetry valid for the in- 
finitely heavy mass of charm quark. Under this symme- 
try, the strong interaction conserves total angular mo- 
mentum of the light quark, j. In the meanwhile, total 
angular momentum of the light-heavy system, J, should 
be still regarded as a good quantum number of the sys- 
tem, even if the heavy-quark symmetry breaks down. In 
this way, the classification of the charmed mesons can be 
explained in terms of the quantum numbers (L,S, J,j), 
where L and S denote the orbital angular momentum 
between the light and heavy quarks and total spin of the 
system, respectively. The doublets with j = 1/2 and 
L = have been observed over the past two decades 
(~ 1975 — 1994) following the discovery of open charm, 
because these states have relatively narrow widths. Re- 
cently the first observations of the scalar charmed mesons 
(J = 1/2, L = 1) have been reported. The very narrow 
Z?jj(2317) was first discovered in the Dfir channel by 
the BABAR Collaboration [lland its existence was con- 
firmed by CLEO 0, BELLEli and FOCUS Collabo- 
rations. Its mass was commonly measured as 2317 MeV, 
which is approximately 160 MeV below the prediction of 
the very successful quark model for the charmed mesons 
0. The BELLE Collaboration [|| has also reported the 
observation of a rather broad scalar meson Dq(2308), and 
the FOCUS Collaboration reported evidence for broad 
structures in both neutral and charged final states that, if 
interpreted as resonances in the J p = + channel, would 
be the L>°(2407) and the £>+(2403) mesons. While the 
mass of the scalar meson, Dq(2308), observed by BELLE 
Collaboration is also bellow the prediction of ref. |5j (ap- 
proximately 100 MeV), the masses of the states observed 
by FOCUS Collaboration are in complete agreement with 
ref. E). 

The spectroscopy of cq and cs pseudoscalar, vector and 
scalar mesons is drawn in Fig. 1, where the theoretical 
predictions of ref. |f| are represented as solid lines for the 
cs and dashed lines for the cq, and the experimental data 
are represented as triangles for the cs and circles for the 
cq. 
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FIG. 1: Spectroscopy of the pseudoscalar, vector and scalar 
charmed mesons. The theoretical predictions of ref. |f| are rep- 
resented as solid lines for the cs and dashed lines for the cq, 
and the experimental data are represented as triangles for the 
cs and circles for the cq. 



Due to its low mass, the structure of the meson 
Dfj(2317) has been extensively debated. It has been 
interpreted as a cs state 0, U 0, 0, [3 , two-meson 
molecular stat e Il3l 

El 

, D — K- mixing [15j , four-quark 
states ES 113 llSj or a mixture between two-meson and 
four-quark states |lflj . The same analyses would also ap- 
ply to the meson £$(2308). 

In the light sector the idea that the light scalar mesons 
(the isoscalars ct(500), /o(980), the isodublet k(800) Hj 
and the isovector an (980)) could be four-quark bound 
states is not new j^, H^]- Indeed, in a four-quark sce- 
nario, the mass degeneracy of / (980) and an (980) is nat- 
ural, the mass hierarchy pattern of the nonet is under- 
standable, and it is easy to explain why a and k are 
broader than / (980) and an(980). The decays a — > irir, 
K — ► Kit and fa, an — * KK are OZI superallowed with- 
out the need of any gluon exchange, while /n — > 7T7t and 
an — > t]tt are OZI allowed as it is mediated by one gluon 
exchange. Since /o(980) and an (980) are very close to the 
KK threshold, the /o(980) is dominated by the 7T7t state 
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and do (980) is governed by the rjn state. Consequently, 
they are narrower than a and k. 

In this work we use the method of QCD sum rules 
(QCDSR) [23| to study the two-point functions of the 
scalar mesons, D sJ (2317), A>(2308) and D (2405) con- 
sidered as four-quark states. In a recent calculation 
the light scalar mesons were considered as 5-wave bound 
states of a diquark-antidiquark pair. As suggested in 
ref. the diquark was taken to be a spin zero colour 
anti-triplet. We extend this prescription to the charm 
sector and, therefore, the corresponding interpolating 
fields containing zero, one and two strange quarks are: 
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eabc^dec{SaCj5Cb)(qdj5Csl 



(1) 



where q represents the quark u or d according to the 
charge of the meson. Since D s j has one s quark, we 
choose the j% current to have the same quantum numbers 
of D s j, which is supposed to be an isoscalar. However, 
since we are working in the SU(2) limit, the isoscalar and 
isovector states are mass degenerate and, therefore, this 
particular choice has no relevance here. 

The QCDSR for the charmed scalar mesons with n 
strange quarks are constructed from the two-point corre- 
lation function 



n(g) 



d 4 x e^<O|TL7 n (a0ii(O)]|O>. (2) 



In the OPE side we work at leading order and con- 
sider condensates up to dimension six. We deal with the 
strange quark as a light one and consider the diagrams 
up to order m s . To keep the charm quark mass finite, 
we use the momentum-space expression for the charm 
quark propagator. We calculate the light quark part of 
the correlation function in the coordinate-space, which is 
then Fourier transformed to the momentum space in D 
dimensions. The resulting light-quark part is combined 
with the charm-quark part before it is dimensionally reg- 
ularized at D = 4. 

We can write the correlation function in the OPE side 
in terms of a dispersion relation: 



n 



OPE 



ds 



(3) 



where the spectral density is given by the imaginary part 
of the correlation function: p(s) = — Im[II (s)]. 

In the phenomenological side, the coupling of the scalar 
meson with n strange quarks, S n , to the scalar current, 
j n , can be parametrized in terms of the meson decay 
constant f s „ as 0: (0\j n \S n ) = 
the phenomenological side of Eq. 



\/2/s„m| n , therefore, 
(0) can be written as 
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where the dots denote higher resonance contributions 
that will be parametrized, as usual, through the intro- 
duction of the continuum threshold parameter s$ |26|. 
After making a Borel transform on both sides, and trans- 
ferring the continuum contribution to the OPE side, the 
sum rule for the scalar meson S n can be written as 



2f s m%e- m sJ M = 



ds e^ M2 p Sn (s) , (5) 



where p Sn (s) = p pert (s) + p m °(s) + p<™>(s) + p^ G ^(s) + 
p mlx (s) + p^ 2 (s) + p< G3 > (s), with 27] 
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(3m 2 -sa), (8) 



which are common to all three resonances and where the 
lower limit of the integrations is given by A = m 2 /s. 
From jo we get: p ms (s) = 0, 
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From ji we get: p ms (s) = 0, 
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Finally from we get 
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FIG. 2: Borel mass dependence of the relative contributions 
of the OPE terms: perturbative (long-dashed line), pertur- 
bative plus quark condensate (dot-dashed line), previous plus 
four-quark condensate (dashed line) , previous plus mixed con- 
densate (solid line). 
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In the numerical analysis of the sum rules, the val- 
ues used for the quark masses and condensates are: 
m s = 0.13 GeV, m c = 1.2 GeV, (qq) = -(0.23) 3 GeV 3 , 
(ss) = 0.8(qq), (qga.Gq) = m%{qq) with m 2 , = 0.8 GeV 2 , 
(g 2 G 2 ) = 0.5 GeV 4 and (g 3 G 3 ) = 0.045 GeV 6 . 
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D, 
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and D 



(2s) 



the scalar charmed 



mesons represented by jo, ji and j2 (in Eq. (JJJ) respec- 
tively. In Fig. 2 we show, as a function of the Borel 
mass, the OPE relative contribution of the: perturba- 
tive (long-dashed line), perturbative plus quark conden- 
sate (dot-dashed line), previous plus four-quark conden- 
sate (dashed line) and previous plus mixed condensate 
(solid line), for the Dq 1s ^ meson. We see that there is 
no good OPE convergence and that the four-quark con- 
densate and the mixed condensate contributions are very 
big, as compared with the perturbative contribution, and 
with opposite signal, in such a way that the final result 
is almost the same as before adding these two contribu- 
tions. One can argue that this is not a good sum rule, 
since there is not a good OPE convergence. We notice, 
however, that this is a common feature of the sum rules 
for currents with more than three quarks. The thre-gluon 
condensate contribution is negligible for all three currents 
and this is why we do not show it in Fig. 2. The same 



behaviour is obtained for D^ s ^ and D\f b> mesons. 

In Fig. 3 we show, as a function of the Borel mass, the 
percentage of the pole and continuum contributions to 

(Is) 

the total contribution for the Z? meson. We see that 
in the Borel window 1.0 GeV 2 < M 2 < 1.4 GeV 2 the 
pole contribution is always bigger than the continuum 
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FIG. 3: Borel mass dependence of the relative contributions 
of the pole (solid line) and continuum (dashed line) contribu- 
tions. 



contribution. Therefore, this is the Borel window that 
we will consider. In order to get rid of the meson decay 
constant and extract the resonance mass, mg,, we first 
take the derivative of Eq. (J5J with respect to 1/M 2 and 
then we divide it by Eq. JSJ to get 



Qds e- s ' M * sp Sn (s) 
£ 2 ds e -s/M* pSn{s) 



(19) 



In Figs. 4 and 5 we show the masses of the D^ s> and 

(Is) 

D resonances, respectively, as a function of the Borel 
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FIG. 4: The D ( ° s) mass as a function of the Borel mass for 
different values of the continuum threshold. Dashed line: 
y/so = 2.6 GeV; solid line: y^So = 2.7 GeV; dotted line: 
= 2.8 GeV. 

mass for different values of the continuum threshold. The 
results for the D^ s ^ resonance is similar to that for the 
D^ s ^ resonance, as shown in ref.|27|. 
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FIG. 5: The D" J mass as a function of the Borel mass for 
different values of the continuum threshold. Dashed line: 
■yfso = 2.6 GeV; solid line: ^/so = 2.7 GeV; dotted line: 
= 2.8 GeV. 

Comparing these two figures we see that the mass of 
is around 100 MeV smaller than the others, since 
the £ ) q 1s ' ) and D^ s ^ resonance masses are basicaly de- 
generated [2?J. While it is natural to expect that the in- 
clusion of a strange quark would increase the resonance 
mass by around the strange quark mass (as was the case 
when one goes from D^ s ^ to Dq 1s '), it is really interest- 
ing to observe that this does not happen when one goes 



from D^ s ^ to D^ s \ In terms of the OPE contributions, 
we can trace this behavior to the fact that the quark 
condensate term is smaller in Dp 2 ^ than in D^ s ^ (due to 
the change from m c (qq) to m c 6), however the inclusion 
of the term proportional to m s m c (which is not present 
in Dp 1 ^), compensates this decrease. 

Considering the variations on the quark masses and on 
the continuum threshold discussed above, in the Borel 
window considered here our results for the ressonance 
masses are given in Table I. 



Table I: Numerical results for the resonance masses 



resonance 


n (0s) 




n ( 2s ) 


mass (GeV) 


2.22 ±0.15 


2.32 ±0.13 


2.31 ±0.14 



Comparing the results in Table I with the reso- 
nance masses given by BABAR, BELLE and FOCUS: 
D+j(2317), Dg(2308) and L>°' + (2405), we see that we 

(Is) 

can identify the four-quark states represented by Dq 
and D^ s) with the BABAR and BELLE resonances re- 
spectively. However, we do not find a four-quark state 
whose mass is compatible with the FOCUS resonances, 
L»o' + (2405). Therefore, we associate the FOCUS reso- 
nances, Dq' + (2405), with a scalar cq state, since its mass 
is completly in agreement with the predictions of the 
quark model in ref. ||. It is also interesting to point 
out that a mass of about 2.4 GeV is also compatible 
with the the QCD sum rule calculation for a cq scalar 
meson [llj . 

One can argue that while a pole approximation is jus- 
tified for the very narrow BABAR resonance, this may 
not be the case for the rather broad BELLE and FO- 
CUS resonances. To check if the width of the resonances 
could modify the pattern observed in the masses of the 
four-quark states, in ref. the phenomenological side 
of the sum rule, in Eq. (JSJ, was modified through the 
introduction of a Breit-Wigner-type resonance form • 
It was shown that the best agreement between the right- 
hand and left-hand sides of the sum rule is obtained for 
ms ~ 2.2 GeV. Therefore our conclusion is that the 
inclusion of the width does not change the value of the 
mass obtained for the resonance. 

Besides de masses, another important point to under- 
stand the nature of the charmed meson states is their 
corresponding decay width. One can ask how, in the 
present approach, it is possible to obtain a extremely 
narrow width for L>+ 7 (2317), while the L>g(2308) state 
remain fairly wide? To compute the decay width of the 
hadronic decay D^ s ^ — > Dfir , for example, one has to 
study the three-point function 

T„(p,p',q) = J d 4 x d 4 y e i - p '- x e iq - y 

x(0\T{j DB (x)jL(y)jt(0)}\0), (20) 
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where p = p' + q, and the currents for the two pseu- 
doscalar mesons in the vertex are 



1 



Hn = -^(«a7/t75"a-d a 7/t75<y, 



J 5 /( 

3d 



(21) 



In the phenomenological side, the three-point function 
in Eq. H20|) is related with the vertex coupling constant, 
, which is related with the decay width through 



(la) 



the relation: 



167rm?, 

oi 



r(D^ ls) - 



(22) 



where A(a, 6, c) = a 2 + b 2 + c 2 — lab — 2ac — 2bc. 

For the light scalar mesons considered as diquark- an- 
tidiquark states, the study of their decay width using the 
QCD sum rule approach was done in ref. 0. In Table 
II we show the results obtained for the different vertices 
studied in ref. |2J], as well as the experimental values. 

Table II: Numerical results for the coupling constants 



vertex 


<?(GeV) 




<T7T + 7r~ 


3.1 ±0.5 


2.6 ±0.2 




3.6 ±0.3 


4.5 ±0.4 


foK+K- 


1.6 ±0.1 




f TY + n~ 


0.47 ±0.05 


1.6 ±0.8 



From Table II we see that, although not exactly in 
between the experimental error bars, the hadronic cou- 
plings determined from the QCD sum rule calculation are 
consistent with existing experimental data. The biggest 
discrepancy is for gf QV + v - and this can be understood 
since the /o — > 7r + 7r~ decay is mediated by one gluon 
exchange and, therefore, probably in this case a s correc- 
tions could play an important role. In the case of the 
decay /o(ao) —> K + K~ , the coupling can not be experi- 
mentally measured due to the lack of phase space. 

In the case of the decay D^ s ^ — > Dfn°, for an isoscalar 

f Is) 

D^"', in the QCD sum rule approach one only gets a 
result different from zero for the coupling constant, if one 
allows a break in the SU (2) symmetry. In this case, the 
coupling is proportional to the difference of the masses 
of the u and d quarks, and the difference of the u and d 
quark condensates. In a preliminary calculation we got 



0.06 GeV, 



(23) 



which gives a decay width T(D^ S> -> D s tt) ~ 8 KeV. It 
is important to notice that, if we have used a isovector 

(Is) 

current for the D state instead of an isoscalar current, 

we would get T(D { ls) -> D s n) ~ 260 MeV. Therefore, it 
seems possible, in this four-quark scenario, to obtain a ex- 
tremely narrow width for D+ J (2317), while the L»[J(2308) 
state remain fairly wide. 



In conclusion, we have presented a QCD sum rule 
study of the charmed scalar mesons considered as 
diquark-antidiquark states. We found that the masses 
of the BABAR, D+j (2317), and BELLE, £$(2308), res- 
onances can be reproduced by the four-quark states 
(cq)(qs) and (cs)(us) respectively. However, the mass 
of the FOCUS resonance, Dq' + (2405) can not be repro- 
duced in the four-quark state picture considered here. 
Therefore, we interpret it as a normal cq state, since its 
mass is in complete agreement with the predictions of 
the quark model in ref. ||. We also obtain a mass of 
~ 2.2 GeV for a four-quark scalar state (cq)(ud) which 
was not yet observed, and that should be also rather 
broad. 
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